Magnetic-field Induced Screening Effect and Collective Excitations 
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We explicitly construct the fermion propagator in a magnetic field background B to take the lowest 
Landau-level approximation. We analyze the energy and momentum dependence in the polarization 
tensor and discuss the collective excitations. We find there appear two branches of collective modes 
in one of two transverse gauge particles; one represents a massive and attenuated gauge particle and 
the other behaves similar to the zero sound at finite density. 
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Systems in the extreme environment with intense fields 
[l[ are attracting theoretical interest recently because ex- 
perimental developments are making it a real possibility 
to observe exotic phenomena in strong fields. In quantum 
electrodynamics (QED) a pair production is expected un- 
der a strong electric field E, which is called the Schwingcr 
mechanism, and the critical strength of E necessary for 
the electron-positron production is given by an order of 
\eE\ ~ ml, i.e. E ~ rv?Je ~ 10 16 y/cm. Such a strong 
field is still beyond the experimentally reachable limit. 
It is known that optimization of spatial and temporal 
profiles of the pulse laser could significantly lower the 
critical E, which enhances a chance further to realize the 
Schwinger pair production in the laboratory Q ■ 

A magnetic field B is also of paramount importance. 
Of course a time-dependent E would be generally ac- 
companied by B. Moreover, the effect caused by B it- 
self is interesting enough; one well-known example is a 
phenomenon called the magnetic catalysis [1| which fa- 
vors a i?-induced chiral condensate. The birefringence 
in a medium of strong B can also be a sizable effect for 
\eB\ ~ m\, which requires B ~ m 2 e /e ~ 10 13 gauss. This 
is again far beyond the limit of experimentally obtain- 
able intensity, while this critical B could be of order of 
the surface strength on the neutron star or magnetar in 
nature. 

One might have thought that it is difficult to attain 
such large B that any new phenomenon can be realized 
in the laboratory. A new possibility was then brought 
from physics of the relativistic heavy-ion collision exper- 
iment which aims to create a new state of matter out 
of quarks and gluons (called the quark-gluon plasma). 
The heavy-ion collision experiments have turned out to 
be successful particularly at RHIC in BNL and more 
data is coming from LHC in CERN. Because heavy 
ions are positively charged and move at almost the speed 
of light, they can generate a gigantic magnetic field if 
their collision geometry is non-central. It was shown that 
the maximum B reached at RHIC energy is as large as 
\eB\ ~ ml ~ 10 18 gauss This is no longer a QED 

energy scale but should be already relevant to physics 
of quantum chromodynamics (QCD) such as the chiral 
magnetic effect, the chiral vortical effect, the chiral sep- 
aration effect, and so on [3, @j. 

Physical properties of QED and QCD matter in an in- 



tense B field are under active investigations by means 
of analytical and numerical studies. In the Monte-Carlo 
simulation an introduction of strong B is feasible and 
various observables have been measured in the latticc- 
QCD calculation Ml as well as model studies on the QCD 
phase transitions [8( and on the color superconducting 
state of dense quark matter Q . In this work we will ex- 
plore the -B-induced screening effect in an analytical ap- 
proach, which is analogous to the screening in the finite- 
temperature field theory. 

First of all let us write down the fermion propaga- 
tor. This part might look technical but explicit calcu- 
lations are useful to clarify how the momentum conser- 
vation holds. Using Ritus' method we can explic- 
itly solve the free Dirac equation. In the following we 
use Landau's gauge fixing, i.e. A° = A 1 = A 3 = and 
A 2 = Bx, which describes B in the third direction. We 
here assume eB > 0, though we can easily relax this with 
minor modification. We define the projection matrix Q 
as P k (x) = diag[/ fc+ (x), fk-(x), fk+(x), fk-(x)] with the 
Landau wave-functions; 



f k+ (x)=Mx 1 -p 2 /eB) (fc = 0,l,...) 
/ fe _(z) =4>k-i{x 1 - P 2 /eB) (k = 1,2,.. 



(1) 



These are the expressions for eB > and fk+ and 
should be swapped if eB < 0. Here 4>k{ x ) is the wave- 
function for the harmonic oscillator with k quanta. 

One can then prove that the background of A 2 = Bx is 
eliminated with this projection matrix and the momen- 
tum is replaced accordingly as p — >• p = (0, — \ / 2eBk,p 3 ). 
Now that we have the complete basis of the solutions of 
the Dirac equation it is straightforward to write down 
the two-point functions in a standard procedure; 



(ip(x)^(y)} 



dp 2 dp 3 ^ e - ic "»( x 



- y °)+ip 2 (z 2 -y 2 )+ip 3 (z 3 -y 3 ) 



(2tt) 2 ^ 2u p 

x P k (x) (w p7 ° - p • 7 + m)P k (y), (2) 

which represents the particle propagation from y to x. 
We note that lo p = \/ p 2 + m 2 is the particle energy. 
In the same way we get (ip(y)ip(x)) for the anti-particle 
propagation from x to y, from which we can get the re- 
tarded, advanced, and Feynman propagators. 
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FIG. 1. One-loop diagram of the polarization tensor. The 
external momentum q is not Landau quantized, while the 
momenta in the loop are quantized and replaced by p and 
q, which are further reduced to p and q in the LLL approxi- 
mation. 



We cannot, however, express those propagators in a 
compact form using the ie-prescription because is not 
commutative with 7 2 . In the lowest Landau- level (LLL) 
approximation that we adopt hereafter, the term involv- 
ing 7 2 is dropped and the Feynman propagator for ex- 
ample simplifies into 



S^(x,y) = (T^(x)^(y)} 
dp dp 2 dp 3 




| e BK* I -&) a +(w I -&) 3 ]. 



Pp,!^ — m + ie 



(3) 



with the spin projection matrix; Pq = diag(l, 0, 1, 0). 
The Landau zero-mode exists for one spin state only and 
Pq is nothing but the projection matrix onto this allowed 
spin state. In the above the index p, refers to only the 0th 
and 3rd components, namely, p-7 = PpTf* 1 = p°j° — p 3 j 3 . 
One might have thought that p 2 /eB on the exponential 
could be dropped in the strong B limit. Wc must keep 
this term, however, because the integration variable p 2 
could take any large number. In fact, as explained be- 
low, the momentum conservation is guaranteed by the 
p 2 -integration with this p 2 /eB term left. 

We are now ready to go into concrete calculation of 
the polarization tensor , which has been discussed re- 
peatedly in various contexts such as the magnetic catal- 
ysis [l2[, the chiral magnetic effect [l3|, and so on. In 
this sense the computation and the result of the polar- 
ization tensor have been already known. Our emphasis 
in this work is thus put in the physical interpretaion de- 
duced from the energy and momentum dependence in the 
polarization tensor. 

The one-loop contribution is, from the di- 
agram in Fig. [1] given as HP 11 ' (a;, y) = 
— (ie) 2 [7 M 5o(x, y)j t/ So(y, x)] . Here we consider the 
polarization for the photon, and if it is for the gluon the 
coupling constant should be g and the trace in color and 
flavor space is necessary. Wc shall work in momentum 
space to calculate, 



iIF"(fc,g) = 



d 4 xd 4 y e iq - x+ik - y mr(x,y). 



(4) 



Here q and k are incoming momenta from the external 
legs attached to x and y, respectively. Usually, without 
external magnetic field, H^ v (x,y) is a function of x — y 
due to the translational invariance, which can be Fourier 
decomposed into exp[— ip ■ (x — y)] and exp[— ip' ■ (y — x)] 
with loop momenta p and p' (see Fig. [I), so that x and y 
integrations yield (2tt) 4 S^ {q - p + p') and (2n) 4 5^ (k - 
p'+p). Then p' is fixed as p—q from the former constraint 
and the latter becomes (27r) 4 <5^ 4 ^(fc + q), meaning that 
the outgoing momentum — k is exactly balanced with the 
incoming momentum q. 

In the presence of the magnetic field, however, the real- 
ization of the momentum conservation is non-trivial be- 
cause the translational invariance is apparently broken 
(not broken in fact) by the inhomogencous vector poten- 
tial A 2 = Bx and H^(x, y) is no longer a function oix—y 
alone. Although the final result is known (see Eq. (117) 
in Ref. [HJ), it would be instructive to make an explicit 
confirmation in the following. 

For B ^ the integrations with respect to a; , y°, 
x 2 , y 2 , x 3 , y 3 are just the standard procedures leading to 
the delta- function constraints which lead to p'° = p° —q°, 
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k 3 + q 3 = 0. Wc can further carry out the x 1 and y 1 
integrations to have 




(q 1 ) 2 + (q 2 ) 2 + 4igV - 2igY 



(5) 



for the a^-integration. Together with the contribution 
from the ^-integration wc find, 



cxp 



4eB 



{q 1 ) 2 + (k 1 ) 2 +2{q 2 ) 2 +2i(q 1 +k 1 )(2p 2 -q 2 ) 



(6) 

Then, since the fermion propagator does not depend on 
p 2 , it is straightforward to perform the p 2 -integration and 
we get / ^ exph^ + fciyj = 27rS{q i + fc i } . ^. 
This is how we can obtain the momentum conservation 
and the overall factor eB/2n that is exactly the Landau 
degeneracy factor. 

After all, in the dimensional rcgularization, the polar- 
ization tensor can be expressed as 

IF" (A;, 9) = {2Tif5^\k+q)c-^ [{ql)2+{q2)2] W{q) (7) 
with 

^<f \ e 2 \eB 



IF"(g) 



r 



2tt 2 



-I{q/m), 



(8) 



where II 1 " = U 2 " = U^ 1 = W l2 = 0, which is concluded 
from Pq^Pq = Pqj 2 Pq = 0. In the final expression we 
have defined a function I(q/m) by 



I{x) 



2/(1 -y)-x~ 



1 



4sin- 1 (x/2) 
xyfi — x 2 



(9) 
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FIG. 2. I{q/m) as a function of q 2 /m 2 . The solid and the 
dashed lines represent the real-part and the imaginary-part, 
respectively. The mass threshold for particle-anti-particle de- 
cay is located at q 2 /m 2 = 4. 



for x < 2, namely, |g| < 2m. This should be expressed as 



I(x) = 1 - 



x\J x 2 — 4 



In 



Vx 2 - 4 



+ 



(10) 



for x > 2, namely, |q| > 2m with g° > assumed. [Be- 
cause we are considering the retarded self-energy, the sign 
of the imaginary part changes if q° < 0.] An imaginary 
part appears in this case because the gauge particle can 
decay into a fermionic particle and an anti-particle above 
the mass threshold. 

We make a plot in Fig. [2] to show the behavior of the 
function I(q/m). We can clearly notice that the mass 
threshold for particle-anti-particlc decay is located at 
\q\ = 2m where the real-part diverges if it is approached 
from the smaller-g side, while the imaginary-part diverges 
if approached from the larger-g side. In the polarization 
tensor in (3+1) dimensions there would be no such sin- 
gularity at the mass threshold because of the phase-space 
factor that compensates for the divergence. 

In the definition @ we can readily take the limit of 
m — > 0, i.e. x — > oo to obtain I(oo) = 1. There- 
fore, in this limit, the polarization simply becomes 
W v {q) = (g^ - q>*q i? /q 2 )e 2 \eB\/2n 2 . This result has 
an interpretation as the (l+l)-dimensional self-energy 
(gi*" — q^q v j q 2 )m 2 multiplied by the Landau-level degen- 
eracy \eB\/2w, where m 2 = e 2 /it is the "photon mass" 



in the massless Schwinger model |14| . Since the massless 
Schwingcr model is exactly solvable and the dimensional 
reduction is complete as we checked above, this result at 
m = is presumably exact. 

Here we shall make a rather technical remark about 
the momentum integration. The naive momentum inte- 
gration for the one-loop polarization in (1+1) dimensions 
is UV finite. Nevertheless, the final results could depend 
on a choice of the rcgularization procedure. To see this, 
let us take an example of LT 00 (g° = 0, \q\ 0) in the 



m = limit. From Eq. (|5J) we already know that the 
correct answer should be e 2 \eB\j2-K 2 in this limit. If we 
perform the p°-integration first [15| . however, we can eas- 
ily prove, 



lim 



dp°2{p°) 2 -p- (p-q) 



2tt 



p 2 (p — q) 2 



(11) 



Then, one may conclude that n 00 (q° = 0, |g| -> 0) = 
which is obviously different from the correct answer in 
the dimensional regularization. One can argue the com- 
putation of n 33 (q° = 0, \q\ — > 0) in the same way to find a 
discrepancy again. Such subtlety originates from the fact 
that the limit of \q\ — > is not commutative to the limit 
of m 2 — > as is noticeable from Fig. [2j In other words 
the physical properties in the region with \q\ < 2m are 
not accessible at all by the knowledge from the massless 
Schwinger model. We pay our attention to this subtlety 
since it is common to process the p°-integration first in 
finite-T calculations. The gauge invariant results ([SJ im- 
plies that the electric-current susceptibility at vanishing 
momenta is zero, while the naive estimates are not [l5j . 
Moreover, interestingly, a statement similar to Eq. (jTTJ) 
is correct for the matter part involving the Fermi-Dirac 
distribution function at finite T. We can explicitly make 
sure that the matter part with m — 0, q° = 0, and 
|q| -> is vanishing just like Eq. (fTTj) . The finite-T 
extension will be reported elsewhere. 

Now we can write the inverse propagator of gauge par- 
ticles as 



-iD" 1 ^ 



(Q) 



g^q 2 



W v {q) 



+ (gauge fixing terms). 



(12) 



We note that we neglect the exponential factor 
e~^B^ q ) +W ) 1 assuming that the transverse external 
momenta are sufficiently smaller than eB, for which the 
LLL approximation can be justified. Let us clarify what 
collective excitations are contained in the above propa- 
gator. For this purpose we should fix the gauge first and 
locate the pole(s) of Dn(uj qi q) as a function of ui q for 
a given q. Instead we can locate the pole position by 
solving the equation, det[D^ 1 (cj«j, q)] = 0. 

Here we adopt the temporal axial gauge A = so that 
the physical degrees of freedom is intuitively manifest. 
Then, D~j^ is a 3 x 3 matrix in the Lorentz index and we 
can show that 



dctLD^, 1 ] = dct 



-q 2 ^ - gV + 5 a 5^ { ^M 2 (q) 



-(q°)V[i 2 -M 2 (q)].. 



(13) 



where M 2 = (e 2 \eB\/2TT 2 )I(q/m). The first factor (q ) 2 
is from the longitudinal gauge field and this unphysical 
contribution is exactly cancelled by the Faddeev-Popov 
determinant associated with the A = gauge. The 
next factor represents one of two transverse gauge fields 
whose dispersion relation is intact. The final factor is 
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[C0 2 q -(q 3 ) 2 ]/m 2 



FIG. 3. Spectral function which is made dimensionless by m 
We choose e 2 \eB\/(2nm 2 ) = 1. 



the most interesting; one of two transverse gauge fields 
is influenced by the screening effect by Ai 2 (q). 

Solving q 2 — M 2 (q) — 0, we can figure out the energy 
dispersion relation ui q (q) of the collective excitation in 
the same way as reading the plasmon dispersion at finite 
temperature. For intuitive understanding we do not show 
the pole position but discuss the spectral function which 
contains the information of both the energy dispersion 
relation and the decay width. Instead of defining the 
spectral function from D^{q) directly, we define p{ui q ) 
for the transverse mode affected by Ai 2 (q) as if it were 
a scalar field, that is, 



lmM 2 (q) 



tt [q 2 - RcM 2 (q)} 2 + [lmM 2 (q)Y 



(14) 



Figure [3] is the numerical result when we set 
e 2 |eB|/(27rm 2 ) = 1. Actually changing e 2 |eB|/27rm 2 
makes no qualitative difference. We can see that there are 
two separate branches below and above the mass thresh- 
old. In what follows we discuss them in order. 

Let us look at the branch above the threshold, i.e. 

\q\ = Jlij 2 — (q 3 ) 2 > 2m. In this case the collective 



excitation (photon) is attenuated because it can decay 
into a particle and an anti-particle. The width is diver- 
gently large near the threshold, which is already clear 
from the dashed line in Fig. [2] As the transverse mo- 
mentum q 2 = (q 1 ) 2 + (q 2 ) 2 increases the width becomes 
narrower and narrower. In the limit of q 2 ^> m 2 , even- 
tually, the dispersion is reduced to that of massive pho- 
ton with its mass given by e 2 \eB\/2ir 2 . This asymptotic 
situation is, as already mentioned, the limit toward the 
masslcss Schwingcr model. 

We can notice that the spectral height grows larger 
as q\ increases. The spectral sum rule should be pre- 
served with the contribution from another branch below 
the mass threshold. As noticed from Fig. [3] this collec- 
tive mode below the threshold lies in the space- like region 



) without attenuation, which looks like a sound 
mode. Such an observation of the sound-like mode has 
been also reported recently in the strong-coupling regime 
in the AdS / CFT analysis [16[ . Since we are working di- 
agramatically and the physical contents are transparent, 
we can give a clear interpretation for this collective exci- 
tation in our perturbative approach; we propose that this 
should be an analogue of the zero sound in the Fermi liq- 
uid theory 

Because of the dimensional reduction, on the one hand, 
physical properties of the system in a strong magnetic 
field looks like (1+1) dimensional. At finite density, on 
the other hand, physical properties are dominated by the 
interactions near the Fermi surface, so that the dynamics 
is effectively (1+1) dimensional too. [Since the momen- 
tum in a certain direction is as large as the chemical 
potential, the momenta in other directions perpendicu- 
lar to the motion are negligible and the system is lo- 
cally (1+1) dimensional if the Fermi momentum is large 
enough.] In condensed matter physics the zero sound 
appears in Fermi liquid even at zero temperature as a 
result of the particle-hole excitation in pseudo (1+1) di- 
mensions where the polarization has a logarithmic singu- 
larity similar to our Fig. [5J In the present case, the mag- 
netic field instead of the Fermi surface causes the pseudo 
(1+1) dimensionality, and the particle-anti-particle ex- 
citation instead of the particle-hole excitation constitute 
sound-like collectivity in the vector channel. 

Usually the zero sound is a longitudinal sound mode as- 
sociated with the density fluctuation which appears only 
when the Landau parameter Fq in the fermionic inter- 
action is repulsive. In relativistic matter the zero sound 
emerges as a gapless branch of the vector meson [l^|. It 
is remarkable that in our case the zero-sound-like mode 
lies in the transverse degrees of freedom, so that it is 
indistinguishable from the massless photon even though 
the microscopic content is totally different. 

In this work we focus on the B-induced screening ef- 
fects on massless gauge particles only. It would be very 
interesting to think about the screening effects on massive 
vector mesons in the same way as in Ref . [l8[ . We would 
then anticipate that a massless collective mode should ex- 
ist in the space- like region as discussed in Ref. [l8| , which 
is the complete counterpart of the zero sound, which we 
name the "magnetic zero sound" here. Recently it was 
pointed out that the p-mcson superconductor could be 
possible in the vacuum with a superstrong B field [l9l ]. 
It is an important question how this electric superconduc- 
tivity should or should not be modified by the screening 
effects and particularly by the existence of the magnetic 
zero sound. 

In summary, in the LLL approximation, we found that 
one of transverse gauge particles is screened by the mag- 
netic field, which leads to two branches of the collective 
excitations. One has a modified dispersion relation with 
a mass and a width whose behavior is reminiscent of the 
plasmon at finite temperature. The other looks like the 
zero sound at finite density, which can be naturally un- 
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derstood from the pscudo (1+1) dimensional nature un- 
der a strong magnetic field. These modes are relevant 
to the situation with \eB\ ^> m 2 ~ \q\ 2 and convey de- 
tailed information of the birefringence in the vacuum (not 
medium) with a magnetic field background. Because a 
huge \eB\ is expected in the heavy-ion collision, we may 
have a chance to probe these collective modes indeed. If 
we replace the overall coupling constant e 2 in .M 2 by g 2 
for quark loops, it would enhance the chance for obser- 
vation and also we can discuss the physics implication 
of our results to the QCD phenomena such as the chiral 



magnetic effect. Using the NJL-typc models we can also 
discuss the realization of the magnetic zero sound in the 
excitation spectra of vector mesons. These applications 
are interesting future problems. 
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